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I. INTRODUCTION 

The problem of atoms and molecules in a magnetic field is both a classic example of time-independent perturbation 
theory and a vexing challenge in the study of neutron stars and white dwarfs. A sufficiently intense magnetic 
field cannot be treated perturbatively. The spectra and properties of neutron-star atmospheres depend crucially on 
magnetic field. Indeed, in the intense magnetic field of a neutron star B > 10^" G the nucleus rather than the field 
acts as a perturbation. The electron is effectively confined to move along the magnetic field lines. 

This work extends classic analytic work on the one-dimensional hydrogen atom to form the basis of a perturba- 
tive treatment of hydrogen in an intense magnetic field. This analytic treatment yields binding energies for B ^ 10^^ 
G whose accuracy rivals that of the recent exhaustive treatment of hydrogen in an magnetic field by Ruder et al. |^ 
with substantially less computational expense. 

We also present a straightforward numerical treatment of the hydrogen atom, the hydrogen molecular ion and the 
helium atom. The electron wavefunction is expanded in a convenient basis, and the Schrodinger equation may be 
solved approximately by diagonalizing a suitably calculated matrix. The effective potential between the electrons and 
between the electrons and the nuclei may be determined analytically, expediting the calculation dramatically. 

II. THE SINGLE ELECTRON PROBLEM 

We begin with the problem of a single electron bound by the combined field of an atomic nucleus and strong external 
magnetic field. The Hamiltonian for the electron is given by 

where we have assumed that the nucleus is infinitely massive, M is the mass of the electron and P = p — e/c A. 

To derive the Schrodinger equation for the electron, we make the replacement p = — zTiV. We take the magnetic 
field to point in the z-direction and choose the gauge where — Bp/2, Ap = = and obtain 

where 1 denotes the spin and spatial coordinates of the electron i. e. ri , cti . The spin portion of the wavefunction 
decouples from the spatial component; therefore, we take the electron spins antialigned with the magnetic field to 
minimize the total energy, i.e. to calculate the ground state. 

For Z = 0, we recover the equation for a free electron in an external magnetic field which is satisfied by a function 
of the form 

V'nmp.(r) =i?„,„(p,(/.)e'^^'^/'^ (3) 

where 
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where a^r = ^JTi/Mloh = y^!c7|e|S Q, and iFi is the confluent hypergeometric function. 

It is convenient to define a critical field where the energy of the Landau ground state ?iw^f/2 equals the characteristic 
energy of hydrogen e^/ao, where the Bohr radius, ao ~ 0.53A. The transition to the intense magnetic field regime 
(IMF) occurs at |§| 



Bi = 2m^c (1)^ « 4.701 x 10^ G. 
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We will express field strengths in terms oi f3 = B/Bj. 

For Z 0, the complete solution may be expanded in a sum of ^pnmp^ since these form a complete set. However, 
for sufficiently strong fields, one can treat the Coulomb potential as a perturbation and use the ground Landau state 
with the appropriate m quantum number as the first approximation to the radial wavefunction; this is known as the 
adiabatic approximation. 

Equivalcntly, the adiabatic approximation assumes that the Coulomb potential does not effectively mix the Landau 
states, i.e. 
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To determine the validity of the adiabatic approximation we calculate this quantity for the first two Landau states 
and TO = 0, 
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where a « 1/137 is the fine structure constant. We find for /3 — 1000 {B = 4.7 x 10^^ G), that the Coulomb potential 
mixes the Landau states of hydrogen by at most 1.4 %. For stronger fields, the mixing is even less important. 
In the adiabatic approximation, we assume that 



(9) 



where Zmu{z) remains to be determined, v counts the number of nodes in the z wavefunction, and we expect the 
axial wavefunctions to be different for different values of the magnetic quantum number to. Wc will use the notation, 
\Qm.v), to designate the eigenstates. 

For ri = 0, the functions Rnm assume a simple form 
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With these assumptions the functions Z^rn{z) satisfy a one-dimensional Schrodinger equation. 
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Performing the integral yields 
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which for large z approaches —Ze^/z. The Schrodinger equation with this potential is not tractable analytically. We 
can take one of two paths. First, the potential may be replaced by a simpler yet still accurate function. Second, we 
attempt to solve the Schrodinger equation numerically. 



III. ANALYTIC SOLUTION 



A. The Approximate Potential 



The potential given in Eq. ^ for arbitrary m may be approximated to within 30 % over the entire domain by the 
much simpler form 
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The double factorial is defined by (—1)!! = 1 and (2?t.+ 1)!! = (2n+ 1)(2?t.— 1)!!. For large m, ^kman asymptotically 
approaches •\/2|m| + Ian, the mean radius of a Landau orbital. 

As we see from Fig. the relative difference between the two expressions is largest near z — kmO-H- For m = 0, 
the difference is greater than 5 % from z = 0.1 to z = 10. We do not expect this approximation to yield eigenvalues 
accurate to better than ~ 10% for wavefunctions peaked in this range. 

We obtain the following eigenvalue equation with the approximated potential. 
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This equation is nearly identical to the Schrodinger equation with a Coulomb potential; therefore, we treat it as a 
Coulomb problem by using the natural units (Bohr radii for length and Rydbergs for energy), 
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where fc„i ^1/2/3 

Again as in the Coulomb problem, we perform the following substitutions 
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This equation may be solved in terms of Whittaker's functions 0. First, we have 
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ZliO = A±M^zs/2m + U) = A±i\^\ + 6n) iFi(l - nZ, 2, |^| + C™)e-^l«l+^"^/' (22) 

where A± are the normahzation constants for ^ > and ^ < respectively. Unless nZ is an integer, these solutions 
tend to infinity as ^ goes to infinity. 

As with the equation for an unmagnetized Coulomb potential, there exists an additional set of solutions. For the 
three-dimensional Coulomb problem, this solution diverges at the origin and is unphysical. However, here we can 
obtain a well behaved solution. By the method of reduction of order, we obtain the alternative solutions, 

Z2{0=A±W,,z,i/2m+U) = A±m+U ii^i(l-nZ,2,|e|+U)e-^l«l+«"^/' / — — 

J [t iFi{l - nZ,2,t)) 

(23) 

These solutions agree with earlier treatments of the one-dimensional hydrogen atom For integer values of nZ, 

the integral in Eq. ^ diverges; therefore, the eigenvalues differ from those of the unmagnetized Coulomb potential. 
Additionally for the unmagnetized Coulomb potential, = and the prefactor of (|^| -I- ^m) is absent. We find that 
in this case, this wavefunction diverges as ^"^ near the origin and only the counterparts of the states given by Eq. |2^ 
are physical. 

The solutions to Eq. ^ will be made of a linear combination of Zi and Z2 ■ For a given magnetic quantum number 
TO, the excitations along the rnagnetic field axis will be denoted by v with u = Q being the ground state. Determining 
the ground eigenvalue of Eq. ^ for a given value of proceeds in reverse. Since the ground state is even, we have 
= A- and Z'{0) — 0. One first selects a value for < nZ < 1. To have the correct behavior as z ^ 00, we 
perform the integral of Eq. H] from |^| + to 00 and calculate Zmo{0 fo'' — 0- 

With the calculated function, one can determine where .^mo(?) ~ ^ ^'^'^ ^^^^ '"^ value of corresponding 
to the eigenvalue nZ. The value of is simply related to the field strength, 

/3=^(-A) (24) 



B. First-order Binding Energies 



As an example we take Z = \ and n — l/-\/l5.58. This corresponds to a bound state (|000)) with an energy of 
15.58 Ry. We find ^0 = 0.141 which yields P = 1000. For /3 = 1000, Ruder et al. |] obtain a binding energy for the 
m = — Q state of 18.60986 Ry. However, it is straightforward to improve upon our estimate of the binding energy 
by treating the small difference between the approximate and effective potential as a perturbation. We obtain 

= (Z„,o|if'|^™o) (25) 

where H' = Vcs — Vapprox- We then obtain the binding energy to first order of 18.48 Ry for (3 — 1000. 

This technique may also be applied to states with m < by using the appropriate value for km in Eq. |l^. For 
example, also for /3 = 1000 and to = — 1 (|0 — 10)), we obtain the zeroth order binding energy of 10.45 Ry and the 
first-order corrected value of 13.71 Ry compared to the result of Ruder et al. ^ of 13.90394 Ry. Since Eq. 15 is 
a better approximation to the effective potential for electrons in the to = state than in to > states we obtain 
eigenvalues to first order within 0.7 % of the fully numerical treatment for (3 > 1000 for these states (compared to 
within 1.4 % for TO = —1 states). 

To calculate the wavefunctions with z/ > 0, we calculate 2^2 (C) for '^■^ > 1 use the first extremum or zero of 
■^2(0 as the value of ^0 for the even and odd solutions respectively. Fig. || depicts .^2(C) for several values of nZ. 
For nZ between k and fc + 1, Zi(x) has k zeros and fc -I- 1 extrema. Therefore, we find that the v > Q states have 
zeroth-order binding energies of fractions of a Rydberg. The calculation of Z2 (■f ) is complicated by the fact that the 
function Z\{^ also has zeros in the range of integration from ^ to cx3 which make Eq. |2^ ill defined. To pass over the 
singularities in the integrand, we integrate the differential equation ^ directly. 

For smaller values of nZ in the range fc to fc -I- 1, the first zeros and extrema approach ^ = 0. Therefore, for larger 
values of /3, the zeroth order eigenvalues of the v > Q spectrum approach the Bohr energies. The energies of the odd 
states approach the Bohr energies from below (i.e. they are more weakly bound), and the even states with the same 
number of nodes are yet more weakly bound |^ . 

Our first-order adiabatic approximation is less accurate for smaller field strengths. For /3 — 100 and to = (|000)), 
we obtain a first-order corrected eigenvalue of 9.348 Ry compared to the numerically derived value of 9.4531 Ry (a 
difference of 1.1 %). However, for fields of i3 > 5 x 10^^ G, the wavefunctions and binding energies derived in this 
section for to < 3 and arbitrary v are sufficiently accurate for all but the most precise analyses. 
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C. Perturbed Wavefunctions 



To obtain first order corrections to the wavefunctions Zmv and second order corrections to the binding energies, 
we fohow the standard techniques for time-independent perturbation theory We must calculate the following 
quantities 

H'^^ = {Z^,\H'\Z^^) (26) 

for a particular value of /?. Since both V^s and V^approx are symmetric about z = 0, iJ^^ is zero for v odd and /i even. 
We obtain 

^(1) - V f27) 

and 

£;(2) - V ^^^"^^^ f28) 

-^"^ - (0) (0) ^^"^1 

For /3 = 1000, the mixing among the v states is on the order of a few percent. The second order corrections to the 
binding energies for the ground {v = 0) state is 10"'^ times the first order correction. For the excited states with v <Q 
the second order correction is less than six percent of the first-order correction; we quote the binding energies to first 
order for the several of the most bound levels of hydrogen for (3 > 1000 in Table || and depict the wavefunctions to 
zeroth order for /3 = 1000 in Fig. ||. 



IV. NUMERICAL SOLUTION 



A. The Basis Set 



We can make substantial progress by carefully selecting a basis to expand the solutions Z^m- Specifically, we choose 

oc 

Zi_,m{z) = 'S^^mkQkiz) (29) 



where 
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Hk{z) are the Hermite polynomials which are orthogonal on the interval —oo to oo with the weighting function 
exp(— z^). The Gk are the solutions to the Schrodinger equation for a harmonic oscillator potential; consequently, 
they provide a complete set for expanding the functions Zi,m{z). 
To obtain the coefficients in the expansion, we calculate the matrix 

Mki = {QklH^igi) (31) 

which is a function of az and the azimuthal state given by m. We calculate this matrix for k,l < N {N = 5 — 50) and 
diagonalizc it. The eigenvalues of this matrix (X^) are E^,n, and the eigenvectors are the coefficients A^mk in Eq- 29. 
Additionally, the functions Z^m{z) and Gkiz) have definite parity; consequently, for even parity solutions to Eq. IJ, 
only the elements of Mki with k and I even need to be calculated. This reduces the size of the matrix from N'^ to 

Because the number of basis functions used is not infinite, we cannot expect the expansion to span the Hilbert space 
of solutions to Eq. Ilj. To estimate the solution, we vary az to minimize the eigenvalue Ai^ corresponding to the bound 



state that we are interested in. By using an expansion of the form Eq. 29. the binding energies and wavefunctions 
may be estimated for excited states along the z-axis. 

Although the functions Qk satisfy a much different equation from Eq. |l^, if sufficiently many Gauss-Hermite 
functions are included, we can obtain highly accurate eigenvalues and eigenvectors. For the ground state (|000)) 
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with the first 31 Gk, we obtain a binding energy of 18.5579 Ry at /3 = 1000, within a factor of 3 x 10"'^ of the result 
of Ruder et al. , 18.60986 Ry. The results are equally accurate for the first excited state (lOOl)); however, states with 
more nodes require additional terms in the expansion to achieve the same accuracy. Fig. |j compares the zeroth-order 
analytic wavefunction with the numerical wavefunction determined by this technique. 

Obtaining an additional few parts per thousand in accuracy can only justify a portion of the additional computation 
required for this numerical technique; however, this technique may be applied to solve the Schrodinger equation for 
potentials more complicated than Eq. nil 



B. The HJ molecule 

Before proceeding to the multi-electron problem, we study the binding energy of the molecule in an intense 
magnetic field. We restrict our attention to the case where the axis of the molecule is aligned with the magnetic field 
direction. This system retains the symmetry under parity of hydrogen, so the numerical technique may be applied 
directly with only two alterations. 

The effective potential is now given by 

Kff,Om,H+(^) = ^off;0m(2 + a) + VcSfirn{z - o) (32) 

and we must vary the internuclear separation 2a to find the minimum binding energy for the entire system (the 
Born-Oppenheimer approximation, e.g. Q). We find the ground state, |000), at (3 — 1000 has a binding energy of 
28.3457 Ry, compared to the LeGuillot & Zinn- Justin ||] result of 28.362 Ry. The internuclear separation is 0.1818ao; 
LeGuillot & Zinn- Justin |^ find 0.181ao. Fig. ^ depicts the wavefunctions of the ground and first excited state |0— 10) 
for H+. 

The accuracy of our analysis is insufhcient to determine if the ungerade state is slightly bound or unbound relative 
to a hydrogen atom plus a proton. However, in the magnetic case, the electron may be excited into the |0m0) states 
whose axial wavefunctions are similar to that of the |000) state. The |0 — 10) state is much less bound at 20.4252 Ry 
than the |000) state (compared to 18.5579 Ry for the H-|-p system). For stronger fields, the |0 — 20) and more excited 
states are bound relative to the H-l-p system. 

Table ^ depicts the numerical results for the ground and first excited state of in an intense magnetic field. The 
ratio of the binding energies of the |000) and |0 — 10) for is approximately equal to the ratio the energies of the 
same states of hydrogen and the same magnetic field strength. This observation provides a quick way to estimate the 
energies of the excited states of from the binding energy of the ground state. 



Table [11 presents results calculated for different values of magnetic field. Our values differ by less than 0.5 % for 
B > 10^^} and by ~ 1% for the two weaker fields considered. We see that the first excited state of becomes 



bound relative to hydrogen atom and a proton at B « 10^^ G. Furthermore, a comparison of Table III with Table 3 
of Lopez et al. |l^ shows that the ungerade state is unbound for B > 10^^ G. 



V. THE MULTIPLE ELECTRON PROBLEM 



A. Approach and Results 

To calculate the atomic structure of multi-electron atoms, we employ a single-configuration Hartree-Fock technique. 
Specifically, we assume that the multi-electron wavefunction is given by a single Slater determinant of one-electron 
wavefunctions. These wavefunctions are varied to minimize the total energy of the system given the constraint that 
each one-electron wavefunction remains normalized. 

This minimization results in the following eigenvalue equations for the individual wavefunctions, 

F(l)V',(l) = e.V'«(l) (33) 

where 1 denotes the spin and spatial coordinates of the electron i.e. ri, (Ti. 
The operator F(l) is the sum of a kinetic and potential energy term 

F{l)^Ho{l) + Vil) (34) 

where the kinetic term is given by the one-particle Schrodinger equation of an electron in the Coulomb field of the 
nucleus. 
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The potential energy consists of a direct and exchange interaction with the other electrons 



where 



ri2 

2 
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Rather than solve the eigenvalue equations directly, we calculate the total energy of the system given a set of wave- 
functions and minimize this energy by varying the parameters of the wavofunctions. 

In a sufHciently strong magnetic field, these equations for the atomic structure become approximately separable in 
cylindrical coordinates. With this in mind, we take the trial wavefunctions to be of the form 



(38) 



Since we are looking for the ground state of these atoms, we take all the electron spins antialigned with the field and 
the radial wavefunction to be given by n = Landau states with each electron occupying a different m state. We 
obtain 



where 
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and ttH = y^hc/\e\B. Wc suppress the spin portion of the wavefunction and use the natural length and energy units 
of the problem an, e^/ciH- 

The total energy of the system is given by 



i; = ^(V.(i)|F(i)|^,(i)). 



(41) 



To expedite the calculation we can integrate over the known wavefunctions in the p and ^ coordinates. Specifically, 
we begin with the integral over (j) in the potential energy terms 
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+ P2 + i^l - -22)^ - 2piP2 C0S((/)1 - 02) 

The expressions for the functions / and g may be simplified by the substitution 4>\ — (j)2 = 2(/3 + 7r/2) and the definition 
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e = 
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resulting in 
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- k"^ sin^ /3 

where F fc) is the complete Legendre elliptic integral of the first kind. The imaginary portion of the integral for 
g must be zero since the Hamiltonian is hermitian (i.e. unitarity). This may be seen by expanding the denominator 
in powers of sin^ (3 and multiplying this series by i sin(2(mi — mj){(3 + tt/2)). The integral of this term is zero. 

Numerical Recipes provides routines to efficiently calculate F (-1, k) = cel(-\/l — k^, 1, 1, 1) (unfortunately, this 
routine is absent from the latest edition). Furthermore, for \mi — mj \ = 1, we can use the same routine to calculate g, 
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g{±l,pi,P2,zi - Z2) = 
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For I mi — rrijl > 1, we must perform the integral numerically. 

Furthermore, we can gain insight on both the functions / and g by expanding them in the limit of large A2 ~ \zi~Zj\ 
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ff(±3, Pi, P2, Zl - Z2) = (27r)^ 

g(±4, pi, p2, Zl - Z2) = (27r)^ 
and in general 
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to leading order in 1/Az. 

In the limit of large Az, the integrals over the radial wavefunctions may evaluated using these expansions. This 
calculation is speeded by the observation that 
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2npdpRo^Ap)RomMp'' = J | ' r + + " + 1 (60) 

y |mi|!|m2|! V 2 J 

which may be proven by using the normahzation condition of the functions Romip) and analyticaUy continuing the 
factorial function with the Gamma function. For Az < 10 we have numerically integrated the functions / and g over 
the various pairs of Landau states. 

After the integration over the radial and angular coordinates, the energy may now be written as expectation values 
of operators acting on the Z(z) wavefunction. Since each electron is assumed to be in a particular Landau m level, 
we can calculate an effective potential energy between the electron and the nucleus by integrating over p, (j). The 
potential is given by Eq. 

The calculational strategy is similar to the single electron case. The quantum numbers m for each electron are 
chosen ahead of time, and the wavefunction Z{z) is expanded as Eq. ^ with each electron having is own variable 
value of az- For each electron j, the matrix 

{M,)m = {Qk\F{i)\gi) (61) 

is calculated. 

The added complication is that the diagonalization of the matrices Mi must proceed iteratively. For the given 
values of a^, the matrices are first calculated assuming that the other electrons (j ^ i) have Ak = 1 for k — Vj. Then 
each electron's matrix is diagonalized and the z/^th eigenvector is used to calculate the interelectron potential for the 
next iteration. The matrices converge after ^5 — 10 iterations. Next, the values of az for each electron are varied to 
minimize the total energy of the configuration. 

For brevity, we discuss the ground state energies and wavefunctions for H2, He and HHe as a function of field 
strength for f} > 1000. For the molecules we again take the molecular axis to be aligned with the magnetic field 
direction. Since we are interested in the ground states of these species we set 1/ = for all the electrons and assign the 
electrons consecutive m quantum numbers beginning with m — 0. Because none of the electrons have axial excitations, 
we are interested in only the most negative eigenvalue of the electron matrices. This eigenvalue is more efficiently 
determined by varying the coefficients in Eq. ^ directly than by diagonalizing the electron matrices iteratively. 

Table IV gives the binding energies of the most tightly bound states of H2, He, HHe and H calculated numerically 
using the variational method. The energies for H are within 1.1 % of the values quoted by Ruder et al. Q for weakest 
field strength common between the two studies. For the stronger fields, the agreement is even closer. For He the 
energies are within 2.5 % of the values of Ruder et al. for the fields that overlap. 

We also computed the binding energies of H2 and H^ and compared the results with the values found by Lai et al. 
We interpolated the Lai et al. results using a cubic spline with ln/3 as the independent variable. The binding 
energies for H2 were within 0.6 — 3% of the Lai et al. results. The agreement for H~ was poorer ranging from 2 — 7%. 
The results for H and H+ (Table H agree to within 0.9% of the Lai et al. values. 

We compare these interpolated results with the results for the species He, HHe and HHe+ to find that the reaction 

H + He ^ HHe (62) 

is exothermic over the range of field strengths considered. However, if there is sufficient hydrogen present, the species 
HHe would quickly be consumed by the exothermic reaction 

H + HHe ^ He + H2 (63) 

for these field strengths, unless 

H + HHe -> H2He (64) 

is sufficiently exothermic (the binding energy and configuration of H2He is beyond the scope of this paper). The 
potential production channels for HHe+, 

He + H2 ^ HHc+ + H" , (65) 

He + HJ ^ HHc+ + H and (66) 

He + H ^ HHc+ + e" (67) 

are endothermic over the range of field strengths considered. We therefore conclude that at least for hydrogen and 
helium, atoms in an intense magnetic field are far more cohesive than adhesive. 
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B. Validity of the Born-Oppenheimer Approximation 



When studying molecules in a intense magnetic field, we have assumed that the nuclear motion and the electronic 
motion decouple, i.e. the Born-Oppenheimer approximation. Schmelcher et al. have examined the validity of 
this approximation in the presence of a strong magnetic field. They performed a pseudoseparation of the equations 
of motion and derive electronic and nuclear equations of motion. Because we have assumed throughout that the 
nuclei are infinitely massive and that the molecules are aligned with the magnetic field, the corrections to the Born- 
Oppenheimer approximation may be neglected. However, the techniques outlined here, specifically the use of one- 
dimensional Coulomb wavefunctions and Gauss-Hermite functions as a convenient and compact basis for the electronic 
wavefunctions of atoms and molecules in intense magnetic fields, can be extended to the more general case where 
these restrictions have been relaxed. 



VI. CONCLUSIONS 



We have developed both an analytic and a convenient numerical technique to accurately calculate the properties 
of simple atoms and molecules in an intense magnetic field. The calculations presented here complement the earlier 
work. We examine two compounds (HHe and HHe+) in addition to the species studied earlier which may form in the 
intensely magnetized outer layers of a neutron star. Additionally, our technique finds both tightly bound and excited 
states efficiently and accurately which is necessary to calculate the radiative transfer of the neutron star atmosphere. 

The techniques presented in this paper complement the recent work in this area. They provide moderately high 
precision with little computational or algebraic expense. Most recent work has focussed on extremely high precision by 
using a Hartree-Fock-like method to reduce the three-dimensional problem to three coupled one-dimensional problems. 
Generally, two of the one-dimensional equations are solved over an functional basis, i.e. Legendre polynomials IT^ ] or 
spherical harmonics [l^ ], and the radial differential equation is solved numerically over a pointlike basis. Fassbinder 
et al. [|l5| and Shertzer et al. [^,0 use a finite element method throughout. 

The spirit of the work presented here is different. We have solved for the wavefunctions using a basis for all 
three coordinates: the Landau wavefunctions in the angular and radial direction and the one-dimensional Coulomb 
wavefunctions or the Gauss-Hermite functions along the axis of the magnetic field. The power of this technique is 
that the basis functions resemble the actual wavefunctions and preserve the symmetries of the potential; consequently 
only a few basis functions (~ 1 — 2 and 20 respectively) are require to reach moderately high precision. 

The work of Kravchenko et al. ^8 19 1 takes an orthogonal approach and achieves very high precision by solving the 
general problem of a hydrogen atom in an arbitrarily strong magnetic field with a double power series in sin 9 and r. 
It remains to be seen whether this simple and accurate technique can be applied to more general problems. 

The properties of the lowest density layers of a neutron star's crust determine the spectral characteristics of the 
radiation expected from the star. One possibility is that linear chains of atoms form in the sur face layers |^-||,|, 
and the atmosphere condenses at finite density. We find that the reactions between hydrogen and helium are unlikely 
to affect the formation of hydrogen or helium chains unless the formation of hydrogen-helium hybrid chains is favored. 

If the envelope is truncated at sufficiently high density, the thermal isolation of the core can be substantially 
reduced [25|. Furthermore, the composition of the outermost layers determines the spectra from the neutron star 
{e.g. p^-2£ I). Without understanding magnetized chemistry in neutron-star atmospheres, is difficult to interpret 
observations of these objects. 



[1] R. Loudon, Am. J. Phys. 27, 649 (1959). 

[2] L. K. Haines and D. H. Roberts, Am. Journ. Phys. 37, 1145 (1969). 

[3] H. Ruder et al, Atoms in Strong Magnetic Fields : Quantum Mechanical Treatment and Applications in Astropfiysics and 

Quantum Ctiaos (Springer- Verlag, New York, 1994). 

[4] L. D. Landau and E. M. Lifshitz, Quantum Mechanics : Non-Relativistic Theory, 3rd ed. (Pergamon, Oxford, 1989). 

[5] V. Canuto and D. C. Kelly, Astr. Sp. Sci. 17, 277 (1972). 

[6] D. Lai, E. E. Salpeter, and S. L. Shapiro, Phys. Rev. A. 45, 4832 (1992). 

[7] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, New York, 1970). 

[8] B. H. Bransden and C. J. Joachain, Introduction to Quantum Mechanics (Longman, Harlow, England, 1989). 

[9] J. C. Le Guillou and J. Zinn-Justin, Ann. Phys. 154, 440 (1984). 



10 



J. C. Lopez, P. Hess, and A. Turbiner, ^stro-ph/9707050 (unpublished) 



W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vettering, Numerical Recipes in C, 1st ed. (Cambridge Univ. 
Press, Cambridge, 1988). 

P. Schmelcher, L. S. Cederbaum, and H.-D. Meyer, Phys. Rev. A 38, 6066 (1988). 
V. Melezhik, Phys. Rev. A 48, 4528 (1993). 

M. D. Jones, G. Ortiz, and D. M. Ceperly, Phys. Rev. A 54, 219 (1996). 

P. Fassbinder et al, Phys. Rev. A 53, 2135 (1996). 

J. Shertzer et al, Phys. Rev. A 39, 3833 (1989). 

J. Shertzer et al, Phys. Rev. A 40, 4777 (1990). 

Y. Kravchenko et al, Phys. Rev. A 54, 287 (1996). 

Y. Kravchenko et al, Phys. Rev. Lett. 77, 619 (1996). 

M. A. Ruderman, in Physics of Dense Matter, edited by C. J. Hansen (Reidel, Dordrecht, Holland, 1974). 
H.-H. Chen, M. A. Ruderman, and P. G. Sutherland, ApJ 191, 473 (1974). 
E. G. Flowers et al, ApJ 215, 291 (1977). 
E. Miiller, A&A 130, 415 (1984). 

D. Neuhauser, S. E. Koonin, and K. Langanke, Phys. Rev. A 36, 4163 (1987). 
L. Hernquist, MNRAS 213, 313 (1985). 

G. G. Pavlov, Y. A. Shibanov, J. Ventura, and V. E. Zavlin AfcA 289, 837 ( 1994). 



G. G. Pavlov, V. E. Zavlin, J. Triimper, and R. Neuhauser, astro-ph/9609097 (unpublished). 
V. E. Zavlin, G. G. Pavlov, and Y. A. Shibanov, A& A 315, 141 (1996). 
M. Rajagopal, R. W. Romani, and M. C. Miller, ApJ 479, 347 (1997). 



11 



TABLE I. The zeroth and first-order binding energies of fiydrogen in an intense magnetic fieid in Rydberg units 
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TABLE II. Tlie binding energy of in an intense magnetic field. Tfie values have been derived numerically and the final 
column gives the numerically derived binding energy of the ground state of H for comparison. 
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TABLE III. The binding energy of in an intense magnetic field for comparison with the results of Lopez et al. 
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TABLE IV. The binding energy of He, HHe"*" and HHe in an intense magnetic field. The number in parenthesis gives the 
number of free parameters in each variational model. The electrons occupy the most tightly bound states, |OmO), e.g. |000), 
|0 — 10) and |0 — 20) for HHc. The values have been derived numerically and the final column gives the numerically derived 
binding energy of the ground state of H for comparison. 
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FIG. 1. The relative difference between the effective potential and the approximated potential. The solid line traces the 
difference for the m = state and the dotted line gives the m = —1 state. 
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FIG. 2. The function Z2{0 for = for nZ = 1/2, 3/2, 5/2. 
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